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^ I Abstract 



We investigate existence and uniqueness of duality solutions for a scalar conservation 
law with a nonlocal interaction kernel. Following [3], a notion of duality solution for such 
a nonlinear system is proposed, for which we do not have uniqueness. Then we prove that 
a natural definition of the flux allows to select a solution for which uniqueness holds. 



1 Introduction 

o 



At a continuous level, many physical or biological systems are modelled thanks to scalar con- 
servation laws. In this note we will focus on a weakly nonlinear system of the kind : 

d t p + d x (a(u)p) = 0, d x u = p, (1.1) 

where a is a given smooth function, a G C°(M). This system is complemented with the initial 
data p(t = 0) = p°. We notice that we can rewrite (II. ip as a single equation since we have 
u = H * p where H is the Heaviside function and we recover the so-called non local aggregation 
equation. This model arises in several applications in physics and biology where a self-consistant 
interaction field u governs the evolution of a density of population p. Then u is defined as 
u = —d x <f) where is the interaction potential. For instance, in the modelling of cell movement 
by chemotaxis, is the concentration of some chemical called chemo-attractant (when a is 
non-increasing) or chemo-repellent (when a is non- decreasing) which drives the dynamics of 
individuals (bacteria). In gas dynamics, this model can be derived thanks to a high- field limit 
from the Vlasov-Poisson-Fokker-Planck system [7], a nonincreasing (resp. nondecreasing) a 
corresponds the to the repulsive (resp. attractive) case. 

From a mathematical viewpoint, it is well-known that in the attractive case, i.e. when a is 
non-increasing, finite time blow-up of regular solutions for such system occurs (see e.g. [T] and 
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references therein). Therefore one has to look for solutions p which are measure- valued in space, 
which generates several difficulties, because the velocity a{u) turns out to be discontinuous, so 
that the product in the divergence term is not well defined, and the corresponding flow has to 
be defined cautiously. A recent approach consists in using techniques from optimal transport, 
see jl]. Another possibility is to define a priori the product. For the Vlasov-Poisson system, 
this has been done in [7J. 

The aim of this note is to interpret (12.21) as a linear conservation equation solved in the 
duality sense [2], the product being defined afterwards, following the strategy introduced in [3] 
for pressureless gases. Therefore we recall in the next section the notion of duality solutions 
and some useful results. In Section 3 we state and prove the main result concerning existence 
and uniqueness of duality solutions of system (II. ip . Section 4 is devoted to some examples of 
applications of this result. 

2 Duality solutions for linear equations 

The notion of duality solutions was introduced in [2] to give a sense to linear conservation 
equations 

d t p + d x (bp) = 0, (2.2) 

when the coefficient b can be discontinuous but satisfies the so-called one-sided Lipschitz (OSL) 
condition 

d x b(t, .) < /3(t) for (3 G L\Q, T) in the distribution sense. (2.3) 

Duality solutions are defined as weak solutions, the test functions being specific Lipschitz solu- 
tions to the backward linear transport equation 

d t p + b(t, x)8 x p = 0, p{T, .) = p T G Lip(R) . (2.4) 



Definition 2.1 1. We say that a Lipschitz solution p to (K41 is a reversible solution if p 

is locally constant on the set 

V e = {(t,x) G [0,T] xR; 3 p e G S, p e {t,x) ± o}. 

2. We say that p G C([0, T];Mi oc (M) - a(Mi oc ,C c )) is a duality solution to / TOP if for 
any < r < T, and any reversible solution p to ft2.4\ ) with compact support in x, the 
function t y f p(t, x)p(t, dx) is constant on [0, r]. 

The most important facts for our purpose concerning duality solutions are gathered in the 
following theorem. 

Theorem 2.2 (Bouchut, James JBj) 

1. Given p° G M.i oc (R), under the assumptions $2.3\) . there exists a unique p G C([0, +00 [, Aii oc (R)), 
duality solution to A2.2\) . such that p(0, .) = p°. 

Moreover, if p° is nonnegative, then p(t, •) is nonnegative for a.e. t> 0. And we have the 
mass conservation \p(t,-)\(M.) = |p°|(R), for a.e. t g]0,T[. 
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2. Backward flow and push-forward: the duality solution satisfies 

Vt E [O,T],V0 E C C (R), [ (j)(x)p(t,dx) = [ (j)(X(t,0,x))p°(dx), (2.5) 

Jr Jr 

where the backward flow X is defined as the unique reversible solution to 
d t X + b(t,x)d x X = m]0,s[xR, X(s,s,x) = x. 

3. There exists a bounded Borel function b, called universal representative of b, such that 
b = b almost everywhere, and for any duality solution p, 

dtp + d x (bp) = in the distributional sense. (2-6) 

Remark 2.3 A similar notion of duality solution for the transport equation is available d t u + 
bd x u = 0, and p is a duality solution of 112. 2\) iff u = J x p is a duality solution to transport 
equation (see J^j). 

We shall need also the following result whose proof can be found in [3] (Theorems 3.1 and 
3.2) 

Theorem 2.4 Let f E C 1 (IR). Let M be an entropy solution to the conservation equation 

d t M + d x f(M) = 0, 
with nondecreasing initial datum M°. Then p := d x M is a duality solution to 

d t p + d x (bp) = 

where we can choose b = f'(M) a.e. Moreover, for all t e]0, +oo[, d x b < 1/t and its universal 
representative b satisfies d x f(M) = bd x M. 

3 Duality solutions for weakly nonlinear equations 

We introduce the following notion of duality solution for the coupled system ( II. ID . inspired by 
the strategy used in [3] for pressureless gases (see also section I4T21 below) . 

Definition 3.1 We say that (p, u) is a duality solution of on ]0, T[ if there exists a bounded 
Borel function b with d x b < a £ L^ oc (0, T) such that 

1. for all < t\ < t 2 < T, d t p + d x {bp) = in the sense of duality on ]t\, ^2 [, 

2. we have d x u = p in the weak sense, 

3. b = a{u) almost everywhere. 
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We underline at once the fact that this definition does not lead to uniqueness as it stands. 
Indeed, assume that a is a non- increasing C 1 function on R and take for initial data p° = 5 X0 , a 
Dirac measure in Xq G R. Looking for a solution as a Dirac mass 5 Xl u), thanks to Remark l2.3l we 
solve the transport equation with coefficient a(H{x — Xi(t))), where H denotes the Heaviside 
function. Then (5 Xl (t), H(x — 2i(i))) is a duality solution of (jl.ljl in the sense of Definition 13.11 
provided that £i(0) = xq and that the admissibility condition a(l) < x[(t) < a(0) holds. Thus 
we have an infinite family of duality solutions. 

Therefore the main result of this note is to explain how a more precise description of the 
product bp in the scalar conservation equation allows to recover uniqueness. It is actually given 
in a very naive way by writing a(u)p = a(u)d x u = d x A(u), where A is an antiderivative of a. 
This choice can be justified in a more rigourous way when the system (II. ip is obtained as the 
hydrodynamic limit of a kinetic system, as it is the case both in [6] and [7|. It turns out that the 
previous formal computation is correct at the kinetic level, so that the flux J e := J £/ E (£) d£, 
where f £ is the distribution function of particles, actually converges to J = d x A(u), which 
defines the flux of the conservation equation. The point now is to justify that this can be used 
to solve the conservation equation in the duality sense. 

Theorem 3.2 Let p° e Mi 0C (R), p° > 0. There exists a unique duality solution (p, u) to the 
non local interaction equation U.l\) in the sense of Definition ^. 1\ which satisfies bp = d x (A(u)) 
where A is an antiderivative of a. 

Moreover, if a is a non-increasing function (i.e. in the attractive case), there exists a flow 
associated to a(u), in other words, there exists a Lip schitz function X such that p(t) = X(t)#p° . 

Proof. Let us denote by u the entropy solution of equation 

8 t u + d x A{u) = 0, (3.7) 

with initial data u° := j x p°(dx). From Theorem 12. 4[ p := d x u is a duality solution of dtp + 
d x (bp) = where we can choose b = a(u) almost everywhere and it is the unique solution 
satisfying bp = d x A(u). Indeed if we have two such nonnegative solutions p\ = d x u\ and 
P2 = d x U2, then u\ and U2 are monotonous solutions of ( 13. 7p . Thus they are entropy solutions 
of this scalar equation and u\ —u^- 

In the attractive case, the entropy solution u is nondecreasing. Therefore p := d x u is 
nonnegative and d x a(u) = a'(u)d x u < since in the attractive case a is non-increasing. Thus 
the velocity field a(u) satisfies the OSL condition (12. 3p and from Theorem 12.21 (ii) there exist a 
backward flow X such that (12. 5p is satisfied. 

In the general case, we can apply the classical Oleinik entropy condition and get that d x b < 
1/t. Then the solution is defined on all ^2 [ f° r < t\ < ti < T and the flow cannot be 
defined up to 0. 

4 Examples 

4.1 Positive chemotaxis 

Equation (II. ip for a non- increasing function a can be obtained from a hydrodynamical limit of 
a kinetic model describing positive chemotaxis (see e.g. [SIE])- Thus from Theorem 13.21 there 
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exists a flow X such that p = X#p°. Let us first come back to the example in subsection 3.1: 
we assume that a is a given non-increasing C 1 function and take p° = 8 X0 . Then we solve the 
Riemann problem 

d t u + d x A(u) = 0, u(t — 0,x) — H(x — x ), 

where A is a concave function. Then the entropy solution is given by u(t,x) = H(x — Xi(t)) 
where the Rankine-Hugoniot condition implies x'^t) = A(l) — A(0). Thus the unique duality 
solution in the sense of Theorem 13.21 is given by (p, u) = (5 Xl ( t ), H(x — Xi(t))) where Xi(t) = 
x + (A(l)-A(0))t. 

On the other hand, if we look for a solution in the form p(t) = 8 X1 (A, then u(t,x) = 
H(x — X\(t)). Integrating equation ( 12. 6ft we get that bp = —dtu = x[(t)5 Xl ^)- By deriving 
in the distribution sense A(u), we get that the definition of the product in Theorem 13.21 bp = 
d x A(u) is satisfied if and only if x[(t) = A(l) — A(0), thus we recover the Rankine-Hugoniot 
condition. Hence the definition of the product allows to select one solution among those found in 
subsection 3.1. It gives more generally the dynamic of aggregates, which are modelled by a sum 
of Dirac masses p° = X] m «<W A similar computation gives the velocity x'^t) = (^4(XT m j) — 
A(Y^~ 1 m j)) / m i- Notice that the velocity of each aggregate is defined by a local equation, 
despite the fact that the initial equation is non local. In the particular case where A is strictly 
concave, aggregates collapse in finite time. 

4.2 High field limit of Vlasov— Poisson Fokker— Planck 

In [7J, the authors prove that solutions to the Vlasov-Poisson-Fokker-Planck system converge 
in the high field limit to solutions of ( 11. II) where a(u) = u in the repulsive case and a(u) = —u 
in the attractive case. To do so, the authors define a weak product pu, which can be proved 
to coincide with the one used here. Applying Theorem 13.21 we can recover the result stated 
in Theorem 2 of [7] : there exists a unique global in time solution of (II. ip in the distribution 
sense such that the product pu = ±u 2 /2. Moreover, in the attractive case, there exists a flow 
X such that p = X#p°, and the dynamics of aggregates is similar to the one of chemotaxis. In 
the general case, the Oleinik entropy condition gives that p < 1/t. Finally, we notice that the 
result of [7] has been extended in two dimensions by Poupaud in [8] by using defect measures 
to define the product of p by u. However, there is no uniqueness of solutions. 

To conclude, we focus on the connection between pressureless gases and the Vlasov-Poisson- 
Fokker-Planck limit which is mentioned in [Tj. The pressureless gases system reads 

dtp + d x {pv) = 0, d t (pv) + d x {pv 2 ) = 0. (4.8) 

Bouchut and James in [3] introduced the notion of duality solution to ( 14. 8ft : 

Definition 4.1 We say that a couple (p,q), p, q £ C([0, T[; Aii oc (M>)), p > 0, is a duality 
solution to ( [^.6] ) if there exists b G L°°(]0, T[xR) and a G Lj oc (]0,T[) satisfying d x b < a in 
]0,T[xR such that 

1. For all < t\ < t 2 < T , we have in the sense of duality on ]t 1) i 2 [xM 

d t p + 8 x (bp) = 0, d t q + d x (bq) = 0; 

2. bp = q. 



5 



The existence result for the Cauchy problem with initial data (p°, q°) strongly exploits the 
relationships between H4.8[) and the conservation law d t u + d x A{u) = 0, where p = d x u, q = 
d x A(u) and A is determined by p and q. Uniqueness follows if A can be defined by the initial 
data, which enforces additional conditions on {p°,q )- 

In the context of Vlasov-Poisson-Fokker-Planck, the function A is given: A(u) = — u 2 /2 
(attractive case) or A{u) = u 2 /2 (repulsive case). Therefore we propose the following variant to 
the results of [3]. 

Theorem 4.2 Let p° G Mi oc (R), p° > and A e C\R). Define q° = d x A(u°), where u° = 
J x p° . Then there exists a duality solution (p°,q°) of the pressureless gases system ^ n ^e 

sense of ^Definition \4-l\ Moreover this solution is the unique duality solution which satisfy the 
relation bp = d x [A(u)}, where u = j x p(dx). 

This theorem is proved in the same way as Theorem 13.21 The solution to the high-field limit 
of the Vlasov-Poisson-Fokker-Planck system obtained in [7] is therefore the unique duality 
solution to ( 14. 8 P given by Theorem 14.21 for the corresponding A. 
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